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ANALYSIS  OF  A PARTICULAR  STOCHASTIC  DIFFERENTIAL  EQUATION 

Ross  D.  Brazee  ^ 


The  stochastic  differential  equation  and  its  analogous  difference  equation  for  a damped, 
oscillatory  system  subject  to  random  impulses  are  discussed.  For  a specific  set  of  conditions 
on  the  coefficients  of  the  differential  equation,  the  autocorrelation  function  and  the  relations 
among  the  coefficients  of  the  differential  and  difference  equations  are  obtained. 


INTRODUCTION 


Bartlett  2 (pp.  144-52,  265-69)  discusses  in  some  detail  the  stochastic  differential  equation 

dX{t)-\- (xX{t)dt-\- ^X{t)dt  = dY{t),  (1) 

and  its  analogous  difference  equation 

Xt+2~Y  aXt+i-\-  bXt=Yt+2>  (2) 

Equation  (1)  may  be  considered  as  representing  a damped,  oscillatory  system  which  is  subjected  to  random 
impulses.  The  notation  dY{t)  represents  the  differential,  in  the  mean  square  sense,  of  the  additive  process 


* Pioneering  Research  Laboratory  on  Physics  of  Fine  Particles,  Agricultural  Engineering  Research  Division,  ARS,  USDA, 
Wooster,  Ohio. 

The  Pioneering  Research  Laboratory  on  Physics  of  Fine  Particles  was  established  in  1962,  at  the  Ohio  Agricultural  Experiment 
Station,  Wooster,  by  the  Agricultural  Research  Service  of  the  United  States  Department  of  Agriculture,  in  recognition  of  the  need  for 
continued  fundamental  research  on  fine  particle  behavior  and,  in  particular,  the  criticality  of  such  studies  to  agriculture. 

In  agriculture,  pesticides  in  both  liquid  droplet  and  sohd  particle  form  are  used  extensively,  but  their  use  is  frequently  hampered 
by  lack  of  optimum  efficiency  and  precision  of  apphcation,  despite  the  best  efforts  of  manufacturers  and  applicators.  This  results  in 
an  excessive  economic  burden  on  agriculture  in  the  maintenance  of  the  best  of  quality  in  agricultural  produce,  as  well  as  an  increased 
chance  of  contaminating  surrounding  areas.  Problems  of  air  pollution  are  also  of  prime  pertinence  and  concern.  It  is  desired  that 
the  results  of  these  investigations  in  fine  particle  behavior  and  its  inherent  subject  matter  will  serve  agriculture  and  its  allied  indus- 
tries, and  other  industries  with  similar  vexations,  in  the  alleviation  of  these  problems. 

^ M.  S.  Bartlett,  An  Introduction  to  Stochastic  Processes.  312  pp.  London.  1956. 
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Y{t)  of  the  accumulated  impulse  effects.  The  impulses  act  to  change  the  velocity  X{t)  discontinuously. 
Bartlett  gives  the  solutions  for  these  equations  as 


ft  pkAt-v)  — p\M-v) 

= + + dY{v) 

j „ Ai  — A2 


(3) 


for  (1)  and 


r-w-2.,  s+l  — 

Z.  = ^ 

s = 0 


/Xl  — )Lt2 


Yr- 


(4) 


for  (2)  Here,  A{u)  and  B[u)  are  determined  from  the  values  of  X{t),  X{t)  at  t—u,  and  A and  B are  obtained 
from  the  values  Xu  and  Xu+i  at  u,  1.  The  values  Xi  and  X2  are  roots  of  the  equation 


\2  + aX  + i8  = 0, 


(5) 


and  /Ti  and  /X2  are  roots  of  the  equation 


+ aEr  b — O, 


(6) 


where  £’r=l  + A is  the  displacement  operator  in  the  calculus  of  finite  differences.  If  |e\|<  1 and  |eSl<  1, 
X{t)  becomes  independent  of  A'(u)  as  u—^  — ^ and  (3)  becomes 


X{t)^ 


p gXi(r  m). 

J ^ 

J —X  ^ 


-u)  — p\2^t-u) 


A.2 


dY{u). 


Similarly,  if  |/xi|  < 1 and  |/u.2|<  1,  then  as  — (4)  becomes 

Y 

J-  r— 


r-w-2;.  s+l  — n^SA-1 


s = 0 


/Lti  — /X2 


(7) 


(8) 


The  above  system  exhibits  a damped  correlogram  when  subjected  to  a time-series  analysis.  The 
coefficients  a and  b may  be  estimated  by  least  squares  methods,  and  if  the  coefficients  a and  (3  are  esti- 
mated from  them,  then  the  theoretical  correlogram  may  be  more  conveniently  calculated.  Bartlett  gives 

the  appropriate  procedure  for  the  case  2=  — In  the  remainder  of  this  paper,  we  shall  discuss  the  case 


THEORY 


Bartlett  ^ indicates  that  the  covariance  function  w{t)  for  (3)  is  given  by 


uKt)  = (Tx^ 


/\ieV  — X.2e^i^\ 

V X1-X2  r 


(9) 


where  cr|  is  the  variance 


ax^  = E{X*it)X{t)},  (10) 

with  X*{t)  representing  the  complex  conjugate  of  the  process  X{t).  Since 

w{T)^(rl  p{t),  (11) 


^ See  page  149  of  footnote  2. 
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(12) 


where  p(r)  is  the  autocorrelation  function  of  X{t),  then 

. XieV  — 

pW= — ttzy: — 


For  the  case  ^ 4/3,  we  obtain  from  (5)  that 


and 

where  y is  defined  as 
We  also  define 


A.i=  — 7 + (y2  — /8)^ 
A.2  = — y — (y^  — /8)^, 
y = a/2. 


\o  = (y2-/3)^  (^o^O). 

/ 

Hence,  upon  insertion  of  (13)  and  (14)  into  (12),  we  find  after  some  rearrangement  that 


which  is  equivalent  to 

p(r)  = sinh  \o|t|  +cosh  XqIt 

\A.o 


pir) 


X. 


+ 


e^o’’  + 


for  a stationary  process. 

We  now  write  the  definition 


p(T)  = Ae  cosh  (Xo|t|— 0). 
Since  we  require  that  p(0)  = l,  then  it  follows  that 

A 

and 

P(t) 


1 

cosh  0, 
cosh  (\q|t|  — 0) 


From  (21)  and  (18),  it  is  possible  to  write 
P(t)  = 


[cosh  Ao|t|  cosh  6 

sinh  Ao|t|  sinh  6 

[ cosh  6 

cosh  6 

cosh  6 


= (cosh  Xo|t|  +tanh  0 sinh  Xo|r|) 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


(19) 

(20) 
(21) 


which  gives 


= e“^|’’|(cosh  Xo|r|  +— sinh  Xo|t|), 
Ao 


d = tanh  ^ 

Ao 


(22) 


Eqs.  (21),  (22),  (15),  and  (16)  afford  a means  of  computing  the  theoretical  correlogram.  We  now  deter- 
mine the  relationships  among  the  coefficients  of  (1)  and  (2)  for  the  case  under  consideration,  viz.,  2=  4/8. 
The  autocorrelation  function  p(r)  satisfies  the  difference  equation 

p(r  + 2h)  + ap(r  + /i)  + bp{r)  = 0,  O).  (23) 

Since  p(0)=l,  we  have  from  (23)  the  requirement 

p{2h)  + ap{h)  -\-b  = 0.  (24) 
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If  (21)  is  inserted  in  (24),  we  find  after  some  manipulation  that 

6““''  cosh  {2koh)  cosh  0 — sinh  (2X.o/i)  sinh  6 

= cosh  {koh)  — b cosh  sinh  (koh)  sinh  6. 

(25) 


Equating  the  coefficients  of  cosh  6 and  sinh  6,  we  obtain  the  pair  of  equations 


cosh  (Xo^)  + ^ e cosh  (2Xo^) 

(26) 

and 

ae~°‘hl2  _ g-a/i  (2koh). 

(27) 

These  equations 

are  easily  solved  to  yield 

a — — 2e~°‘^^^  cosh  (koh) 

(28) 

and 

b = e~^^. 

(29) 

Eqs.  (28)  and  (29)  with  (16)  permit  estimation  of  a and  /3  from  the  least  squares  estimates  of  a and  b. 
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